February 22, 2008 , , 

Prom particles to spins: Eulerian formulation of supercooled liquids and glasses 
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The dynamics of supercooled liquid and glassy systems are usually studied within the Lagrangian 
representation, in which the positions and velocities of distinguishable interacting particles are fol- 
lowed. Within this representation, however, it is difficult to define measures of spatial heterogeneities 
in the dynamics, as particles move in and out of any one given region within long enough times. 
It is also non-transparent how to make connections between the structural glass and the spin glass 
problems within the Lagrangian formulation. We propose an Eulerian formulation of supercooled 
liquids and glasses that allows for a simple connection between particle and spin systems, and that 
permits the study of dynamical heterogeneities within a fixed frame of reference similar to the one 
used for spin glasses. We apply this framework to the study of the dynamics of colloidal particle 
suspensions for packing fractions corresponding to the supercooled and glassy regimes, which are 
probed via confocal microscopy. 
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INTRODUCTION 

The phenomenology of structural and spin glasses has 
much in common: no static long-range order, aging relax- 
ation, heterogeneous dynamics, etc While a precise 
and unambiguous connection between these two prob- 
lems still lacks, the possibility that such relation exists 
dates back to the work by Kirkpatrick, Thirumalai and 
Wolynes H 0, H, who proposed a connection between 
structural glasses and the p-spin disordered model. More 
recently, Tarzia and Moore [5[ have paralleled the phe- 
nomenology of structural glasses to that of an Edwards- 
Anderson model in a uniform magnetic field. One of the 
main hurdles in making a direct real space connection 
between these two problems is that spin glass models are 
defined on a lattice, while the particles comprising struc- 
tural glasses are itinerant. 

Supercooled liquids and glasses are usually described 
within the Lagrangian formulation, where one tracks the 
position of individual particles as a function of time. Nat- 
ural quantities computed within this frame of reference 
are mean-square displacement and self-diffusion of these 
particles. Heterogeneous dynamics can be probed, for 
example, by studying quantities such as mobility within 
prescribed boxes; however, such fixed regions serve this 
purpose just for a certain time, as particles move in and 
out of these boxes if one waits for long enough. In con- 
trast, studying local dynamics in a spin glass presents no 
such complication, as spins remain fixed to their sites at 
all times, and all that changes is the spin orientation as 
function of time. Therefore, if one is to construct a sim- 
ple description of particle systems that could actually be 
used in analyzing real experimental data from the point 



of view of a spin glass, one must abandon the Lagrangian 
formulation. 

We propose an Eulerian analysis of the dynamics of 
interacting particle systems. The proposal consists in 
transforming the data of numerical simulations or con- 
focal microscopy experiments, usually presented in the 
Lagrangian representation as time-dependentpositions 
and velocities of distinguishable particles [H, 0, la, Q , into 
time-dependent occupation numbers of finite volume pix- 
els at fixed positions in space. In this way, we divide 
the simulation or experimental fixed volume box contain- 
ing the particles into pixels located at sites labeled by 
i = 1,...,N, and assign a 'spin', Sj = ±1, depending 
on whether a 'piece' of the particle falls within the pixel 
or not. The data treatment here proposed allows one to 
make close contact with spin (disordered or constrained) 
models. 



METHOD 

We divide the experimental box in cubic pixels of linear 
size a = R/q where R is the radius of the particles and 
q is a parameter, typically with q > 1. The number of 
pixels is then N = V/a d with V the total volume of the 
experimental box (we focus throughout the paper on the 
d = 3 case relevant to the experiments analyzed below). 
The simplest definition of the spin variable is such that 
Si = 1 whenever a particle (independently of which one it 
is) overlaps the i-th pixel, and Si = — 1 otherwise. With 
such a definition, though, the magnetization density is 
non-zero, m = N^ 1 5^ i=1 Sj ^ 0, at a generic volume 
fraction tf>. To work at zero magnetization density and 
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make closer contact with usual spin (glass) problems, we 
shrink the particle size to an effective radius R e g such 
that the covered volume is 50% [lol ]. 

An efficient algorithm that maps particle positions into 
spin variables works as follows. First, one constructs the 
grid of pixels and sets all spins to Sj = — 1 for all i. Next, 
one reads the particle centers from the data file and sets 
to +1 the spin variables of the pixels around the center of 
each particle. One thus avoids having to go over all sites 
in the lattice and to compute distances between particle 
centers. This procedure is repeated at each time step. 

The data set has now been transformed into spin val- 
ues and all interesting correlation functions in the spin 
realization inform us about the dynamics of the particle 
system. The spin variable is naturally related to an oc- 
cupation number, by rij = (sj + l)/2 = 0,1, and then 
to a density. We stress here that, in this construction, 
these densities are not coarse grained quantities built by 
looking at distances larger than the particle size, but in- 
stead the other way around, by looking at distances of the 
order and below the particle size. Within this construc- 
tion, the parallel with the spin-glass problem is also clear: 
a short-ranged equal-time spin-spin correlation function 
corresponds to a short-ranged particle density order, etc. 

ANALYSIS 

We apply this framework to experimental data from 
colloidal suspensions, both in the supercooled liquid 
regime 0, 0] and in the dense glassy phase [§]. The 
suspensions are of colloidal poly-methylmethacrylate 
(PMMA) with radius R = 1.18 fim (and a polydispersity 
of ~ 5%), suspended in a mixture of decalin and either 
cycloheptylbromide (for the samples with (fi < (fi g 0.58) 
or cyclohexylbromide (for the sample with <fi > (fig). 
These solvent mixtures match the index of refraction of 
the particles to aid in visualization. Furthermore, the 
solvent mixtures also match the particle density, so that 
sedimentation does not occur during the experiments. In 
these solvents, the particles are slightly charged, modi- 
fying their pair correlation function somewhat from that 
of hard spheres, although they still undergo the glass 
transition at <fi g ~ 0.58. The particles in dilute samples 
diffuse their own diameter in 11 s, although in these con- 
centrated samples their motion is much slower Q. All 
samples are stirred prior to data acquisition. The two 
samples with (fi < (fi g are stirred to break up any crystals, 
and data acquisition is started after transient flows have 
diminished (~ 30 min). For the sample with <fi > <fi g , 
no crystals are present prior to stirring; instead, stirring 
initiates aging in the sample, and data acquisition begins 
immediately after the stirring is ended, setting the initial 
time t w = [|. 

Data is obtained via confocal microscopy [ll|, which 
is used to rapidly obtain a three-dimensional image of 
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FIG. 1: (Color online) One-time two-spin correlation c(r, t) as 
a function of distance r. Panel (a): evolution with time of the 
glass (<f) ~ 0.62) pair-correlation function. All curves are very 
similar, apart from the last one that does not decay at long 
distances. See the text for a discussion on this fact. Panel 
(b): the time-averaged quantity c(r) for the three packing 
fractions, (fi ~ 0.52, 0.56, and 0.62. 

dimensions approximately 60 x 60 x 12 /im 3 . Within each 
image, particle positions are obtained with an accuracy 
of 30 nm in x and ?/, and 50 nm in z (along the optical 
axis of the microscope). For other experimental details, 
see Refs. 

We choose q = 5 and thus a = R/q ~ 1.18 /zm/5 ~ 
0.236 fim. The 3d positions of the particles were recorded 
every 18 s for the supercooled data sets at <fi ~ 0.52 and 
(fi ~ 0.56, and 20 s for the glassy data set at (fi ~ 0.62. For 
the effective radii we find: R D g = 1.17 /xm at <fi ~ 0.52, 
i? c ff = 1.11 /im at (fi ~ 0.56 and i? e ff = 1-10 /jm at 
0-0.62. 

We now show how to characterize the dynamics of the 
colloidal system of particles using solely the mapped spin 
variables. We start by defining two-spin correlations 

C 2 (r;t,t w ) = ^ Yl . (!) 

i,j;\ri —Fj | — r 

which can be used to determine both equal-time spatial 
correlations and same-site two-time correlations. 
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FIG. 2: (Color online) Two-time two-spin lo- 

cal correlation C{t,t w ) as a function of time- 
delay t — tw in the supercooled liquids with cj> ~ 
0.52, 0.56. The decay at several waiting-times (t w = 
180 s, 306 s, 414 s, 558 s, 756 s, 1008 s, 1350 s, 1800 s) 
is shown with data points connected with thin lines; the 
average of these sets is shown with thick lines, along with 
exponential and stretched exponential fits to the averaged 
data. 



In Fig. [I] we present the equal-time correlation function 
between two spins at a distance r, c(r, i) = Ca(r;t, t), 
which is analogous to the pair correlation function con- 
volved with a square hat function of width i? c ff. Due to 
the finite size of the sample one expects time-dependent 
fluctuations. In the supercooled liquid regime these are 
present but no systematic trend is visible (not shown). 
The time-dependence in the glass is shown in panel 
(a) where the pair correlation function as a function 
of r, computed at equally spaced times, is displayed. 
The curves show no systematic time-dependence until 
t ~ 4000 s. A clear departure is seen at later times when 
the pair correlation no longer decays to zero. The satu- 
ration at long-distances is not related to crystallization 
since Cianci et al found no increase in crystalline order 
as the sample aged [8j . We do not know the exact reason 
for this saturation. In what follows we just analyze glass 
data for times that are shorter than t ~ 4000 s. 

The time averages of the equal-time correlation func- 
tion, c(r) = fc^ 1 Ylk=i c ( r i*fc)i are shown in panel (b) of 
Fig. Q] for the three packing fractions. These are calcu- 
lated using k m = 10 times equally spaced over an interval 
of approximately 6300 s in the supercooled liquid and 8 
times before 4000s in the glass. Notice that the peaks 
move slightly to lower values of r for increasing values of 
<f>, but there in no qualitative difference in this one-time 
quantity for these three packing ratios. These curves are 
essentially the same as the ones shown in Q for the su- 
percooled liquid and in [8j for the glass, computed using 
directly the particle positions. 

We now turn to two-time quantities, starting from 
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FIG. 3: (Color online) Two-time two-spin lo- 

cal correlation C(t,t w ) in the glass. Panel (a): 
data for the decay at several waiting-times (t w = 
180 s, 306 s, 414 s, 558 s, 756 s, 1008 s, 1350 s, 1800 s) 
are shown with thin lines, plotted as function of the time 
delay t — tw', the smoothed decay is highlighted with thick 
lines. Panel (b): scaled data using the simple aging form 
C(t,t w ) ~ f(t/t w ); the solid (black) line is the power law 
/(*)~z- - 36 . 



the global equal-space two-spin correlation C(t,t w ) = 
C2(r = 0;t,t w ). In Fig. [2] we present its decay, as a 
function of t — t w for the supercooled liquid regime. The 
group of curves that fall below are for (f> ~ 0.52. The 
curves drawn with thin lines represent data for several 
waiting-times and, within the numerical error, they have 
the same decay, proving that the dynamics are stationary. 
The thick (red) line is the average over all waiting-times. 
The thick (blue) curve lying above is the averaged data 
for <f> ~ 0.56. The spreading for different waiting-times 
(not shown) is similar to the one for </> ~ 0.52. In both 
cases C decays from 1 to 0.8 in less than 18 s (the mini- 
mum time step for which data is recorded), due to Brow- 
nian motion of the particles within their cages The 
dotted black lines are exponential fits, f(x) = ae~ x / b , 
to the decay for t — t w > 90 s that have been trans- 
lated to make the curve visible. The solid black line is 
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FIG. 4: (Colour online) The four point correlation Ga in the 
supercooled liquid for several time delays given in the key. 
Panel (a) cj> ~ 0.52 and panel (b) <f> ~ 0.56. 



TABLE I: Fitting parameters for global two-time correlation 
decay in the supercooled liquid. The errors (not quoted) are, 
at most, 3%. 
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a 


b(s) 


c 


d e (s) 


0.52 
0.56 


0.75 
0.80 


9950 
23000 


0.75 
0.85 


0.85 10650 
0.65 43400 



a fit of the data for <f> ~ 0.56 to a stretched exponential, 



ce 



-(x/e) a 



fit to the averaged curves. In TableUwe 



give the values of the fitting parameters for both densities 
although in Fig. [2] we show only the stretched exponential 
for the higher packing fraction. 

Figure [3] shows the two-time correlation function in the 
aging regime. Panel (a) displays the relaxation after sev- 
eral waiting-times. As in the supercooled liquid regime, 
the correlation decays from 1 to 0.7 rapidly and then fur- 
ther decays to zero in a much slower manner. In a dou- 
ble logarithmic scale the separation between the station- 
ary (C ~ 0.7) and aging (C ~ 0.7) regimes is seen as a 
plateau at the Edwards- Anderson value q eq ~ 0.7. Panel 
(b) demonstrates that the aging data can be satisfactorily 
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FIG. 5: (Colour online) The four point correlation C4 in 
the supercooled liquid and the glass. Panel (a): compar- 
ison between three packing fractions, <f) ~ 0.52, 0.56 at 
t — t w = 1800 s and the glass <f> ~ 0.62 at t — t w — 2000 s and 
t w = 2000 s. Panel (b): C4 for the glass at several waiting- 
times and t — t w — 1200 s fixed. 



scaled using the 'simple' aging form C(t,t w ) ~ f(t/t w ) 
with f(x) ~ x~ ' 35 for waiting-times that are longer than 
t w ~ 1200 s. However, the range of variation of both axis 
is smaller than a decade and it is hard to give a concrete 
conclusion on 'simple' aging in this sense. Still, it is in- 
teresting to note that this behavior is remarkably similar 
to the one found with Monte Carlo simulations of the 3d 
Edwards- Anderson (EA) spin-glass [121] . 



A two-time dependent correlation length [13j, ll4j can 
be extracted from the spatial decay of a two space points 
and two-times correlation function: 
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i(t)si(t w )sj{t)sj(t w ) ,(2) 



or a variation in which we extract the square of the two- 
time local correlation C(t,t w ) that is the expected large 
distance limit of eq. (J2J): 
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FIG. 6: (Colour online) The time dependence of the corre- 
lation length normalized by the cut-off distance r ma x (= 30 
lattice units) in the supercooled liquid and glass, computed 
using eq. Q . 



These definitions are simple extension of the ones used in 



the analysis of the stationary supercooled liquid 

Figure |4] displays the four-point correlation function 
for several time-delays in the supercooled liquid phase. 
In Fig. [5] we show the space-dependence of the four-point 
correlation function C4 for a fixed time-delay, t — t w — 
1800 s for the supercooled liquid and t — t w = 2000 s for 
the glass. The supercooled liquid curves have been aver- 
aged over the waiting-time taking advantage of station- 
arity. The glassy curve has been smoothed by averaging 
over two time-windows of length r = 200 s centered at 
t w and t. 

The correlation length can be evaluated using 
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£2 _ JO 



dr r 2 Ci{r; ,t, t w ) 
Jo" 10 * dr d(r\t,t w ) 



(4) 



in the limit r max — > 00. This analysis, applied to the 
data in Fig. [5l yields a correlation length of the order 
of £ ~ 4 — 6R. This can be confirmed by simple visual 
inspection since all curves decay close to zero at distances 
r ~ 8 — 10 (im~ 4 — 6R. These values are of the same 
order as the ones found in previous studies 0]. 

In order to capture the time-dependence of £ we use, 
instead, a finite value of 

The reason is that we do not have enough precision at 
t > r max to disentangle the curves measured at differ- 
ent times. We thus obtain slightly shorter values of £ 
that have, though, a systematic temporal trend. Fig- 
ure [H shows these results. The supercooled liquid curves 
follow the expected trend: the sample with a higher pack- 
ing fraction (</> ~ 0.56 with thick green line) has a longer 
correlation length than the one with the lower packing 
fraction (</> ~ 0.52 with thick red line). In both cases 
the length smoothly increases in time. We then com- 
pare these results to the measurements in the glass at 



different waiting-times (thin lines in the same figure). 
All curves grow as a function of time-difference. At 
short time-differences the curves with shorter waiting- 
time have a longer correlation length while the trend re- 
verses at longer time-differences. In the glass the growth 
with time-difference is faster than in the supercooled liq- 
uid and one expects the longer waiting-time curves to 
go beyond the supercooled liquid ones at longer time- 
differences (not reached in the experiment). The reason 
why the correlation lengths in the glass are shorter than 
the ones in the supercooled liquid at the available times 
is that the glass is still far out of equilibrium and corre- 
lations have not propagated far in the sample yet. 

The two-time dependence in the glassy regime is sim- 
ilar to the one found in the 3d EA spin-glass [13], the 
Lennard- Jones mixture [17 1, and the 3d random field 
Ising coarsening system [18 . In the latter case the ori- 
gin of the two-time dependence of the growing length £ 
can be traced back to the one-time dependence of the av- 
eraged radius of the growing domains of two competing 
equilibrium states. In the structural and spin glass cases, 
the two-time dependence of £ does not have such a clear 
simple origin, and it is less well understood. 

We now turn to the study of local correlations, which 
are probes of local heterogeneities in the dynamics. The 
particles in the colloidal system do not displace all at 
the same rate: some regions can reconfigure much faster 
than others, for the same elapsed time between frames. 
A broad distribution characterising these heterogeneities 
can be captured, in the mapped spin system, by using a 
local two-time spin-spin correlation averaged over a cell 
of size V r = £ 3 centered at site r: 



Ct(r;t,t w ) = — ^2 

r i=l£V r 



(5) 



Whenever the cell size P is much larger than the dynami- 
cal correlation length £, the local correlations just reflect 
the global value Ci{r — 0;t,t w ). Instead, whenever the 
coarse graining box I is smaller than £, the local val- 
ues are non-uniform. This fact is captured by a broad 
probability distribution function (PDF) of the local cor- 
relations P(Ct) at fixed times t and t w . 



A simple scaling hypothesis discussed in [12], [14| im- 
plies 



P(Ce 



ft P 



P(Ci;C,t/S,e/L^0,4>) (6) 



with C and £ the values of the global correlation and 
correlation length at the measuring times t and t w , and 
L the size the of the sample that is much longer than the 
coarse graining length. We kept an explicit dependence 
on the control parameter in the system, that is to say, the 
packing fraction 4>. This form is obtained by exchanging 
the dependence on times t and t w by a dependence on 
the two-time dependence quantities C and £ (exploiting 
their monotonicity properties) and then assuming that 
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FIG. 7: (Colour online) The PDFs for local two-time spin-spin 
correlations for (a) supercooled liquid at <f> ~ 0.52, (b) super- 
cooled liquid at (j> ~ 0.56, and (c) the glass at 4> ~ 0.62. The 
distributions signal heterogeneous dynamics within regions of 
linear size I = 2.2R ~ 2.60 /im. Supercooled systems show 
data collapse at fixed time differences t — t m , reflecting time- 
translation invariance, while the glassy sample shows data 
collapse at fixed t/t w , reflecting approximate simple aging 
behavior. 



the three lengths L, £ and £ can only appear through 
the ratios t/£, and t/L. Since the time-variation of the 
correlation length is very slow, as a first approximation 
one can neglect the scaling variable and simply check 
the scaling form P(Cf,C,<f>) 13}. In Fig. [7|we test this 
scaling form using a coarse graining box I — 2.2R ~ 
2.60 /im. The PDFs are shown for different values of the 
waiting time for the two supercooled systems at <fi ~ 0.52 
(a) and <f> ~ 0.56 (b), and for the glass with <f> ~ 0.62 
(c). The PDFs collapse for fixed time difference At = 
t — t w in the case of the supercooled samples, and for 
fixed ratio t/t w in the case of the glass, reflecting that 
time-translation invariance is manifest in the supercooled 
liquid regime, but broken in the glass, which ages. (In 
the case of Fig. [TJb) the bad collapse for t w = 900 s may 
be attributed to lack of equilibration at this high packing 
fraction: fluctuations may be more sensible than average 
values in detecting a remanent time-variation. Another 
explanation would be that the sample we are using is 
too much heterogeneous and not fully representative of 
equilibration at this packing fraction.) As expected the 
PDFs get wider for longer time-differences or larger value 
of t/t w . 

Remarkably, once the PDFs have been scaled, the scal- 
ing function of the sample with low packing fraction is 
very similar to the one of the glass. It is worth noting here 
that the average correlations in the loose supercooled liq- 
uid and the glass are very similar: (C) ~ 0.5. The aver- 
aged two-time correlation in the dense supercooled liquid 
(4> — 0.56) remains, during the available time-window, 
too high to be compared with the other two cases. The 
similarity between the PDFs for <f> ~ 0.52 and 0.62 sug- 
gests that a 'universal' PDF connecting the fluctuations 
for different packing fractions through a proper rescaling 
of times might exist. We plan to explore this hypothesis 
using molecular dynamics of Lennard- Jones mixtures. 



CONCLUSIONS 

In short, we introduced a simple method to translate 
particle data position into fixed frame spin variables. We 
computed correlation functions and extracted a correla- 
tion length from confocal microscopy data of supercooled 
and glassy samples and we found remarkably similar re- 
sults to the ones obtained with numerical simulations of 
spin models. 

One advantage of our analysis method is that it does 
not rely on tracking individual particles over time. In- 
stead, it acts on the particle locations at each time, with- 
out regard to their identity. Thus, this method will be 
useful in situations where particles move too rapidly to 
be tracked, but which are still able to be visualized by 
microscopy in quickly obtained images. 
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